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On the Interpretation of the I — v Features in the Milky Way Galaxy 
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Abstract 

We model the gas dynamics of barred galaxies using a three-dimensional, high-resolution, 7V- 
body+hydrodynamical simulation and apply it to the Milky Way in an attempt to reproduce both the 
large-scale structure and the clumpy morphology observed in Galactic H I and CO I — v diagrams. Owing 
to including the multi-phase interstellar medium, self-gravity, star-formation and supernovae feedback, the 
clumpy morphology, as well as the large-scale features, in observed l — v diagrams are naturally reproduced. 
We identify in our I — v diagrams with a number of not only large-scale peculiar features such as the '3-kpc 
arm', '135-km s _1 arm' and 'Connecting arm' but also clumpy features such as 'Bania clumps', and then 
link these features in a face-on view of our model. We give suggestions on the real structure of the Milky 
Way and on the fate of gas clumps in the central region. 

Key words: Galaxy: disk — Galaxy: kinematics and dynamics — galaxies: ISM — galaxies: spiral 
— method: numerical 



1. Introduction 

In order to understand the structure and dynamics of 
the Milky Way galaxy, the line-of-sight velocity of the 
interstellar medium (ISM) (e.g., Hi, CO) is often used. 
This is used to determine the radial mass distribution 
(Sofue & Rubin 2001) and the spatial distribution (Oort 
et al. 1958; Nakanishi & Sofue 2003; Nakanishi & Sofue 
2006). These studies had to adopt the fundamental as- 
sumption here is that the ISM rotates axisymmetrically in 
a purely circular fashion on the Galactic plane. However, 
longitude- velocity (I — v) diagrams of H I (Hartmann & 
Burton 1997) and CO (Dame et al. 2001) in the in- 
ner Milky Way galaxy are inconsistent with the circular 
motions 1 . Moreover, the central stellar bar and spiral 
arms result in a non-axisymmetric gravitational poten- 
tial. Thereby, there is a large difficulty in reconstruct- 
ing the true mass distribution and spatial distribution of 
the ISM from the line-of-sight velocities ( Koda & Wada 
2002; Gomez 2006; Pohl et al. 2008; Baba et al. 2009 here- 
after Paper I). 

The peculiar features of the I — v diagrams of H I and 
CO in the Milky Way galaxy are widely accepted to be as 
non-circular flows under the influence of the stellar bar 
and spiral arms (see chapter 9 in Binney & Merrifield 
1998). Binney et al. (1991) modeled ISM motion in terms 
of closed stellar orbits (x\ and X2 orbits), and interpreted 
the peculiar features seen in I — v diagrams of the Galactic 



See their Figure 1 and section 2 in Rodriguez-Fernandez & 
Combes (2008) for a recent summary of outstanding peculiar 
features in the l — v diagram of the ISM in the Milky Way galaxy. 



central region (\l\< 10°) as non-circular motions due to the 
Galactic bar. Jenkins & Binney (1994) tried to reproduce 
the I — v diagram of the central molecular zone (CMZ) by 
numerical simulations using the so-called sticky-particle 
method. Many hydrodynamical simulations of the ISM 
have also been performed in order to investigate the ori- 
gin of the H I and CO I — v diagrams of the Galactic disk 
(Wada et al. 1994; Englmaier & Gerhard 1999; Weiner & 
Sellwood 1999; Fux 1999; Bissantz et al. 2003; Rodriguez- 
Fernandez & Combes 2008). These previous numerical 
studies showed that the large-scale features of the I — v 
diagrams largely depend on the mass model of the bulge, 
stellar bar, and spiral arms along with the pattern speed 
and the location of the observer. 

Englmaier & Gerhard (1999) performed hydrodynamic 
simulations of gas flows in the gravitational potential of 
the near-infrared luminosity distribution of the Milky Way 
(Binney et al. 1997). Their best fit models qualitatively 
reproduced number of observed features in the l — v di- 
agrams. Bissantz et al. (2003) extended the model of 
Englmaier & Gerhard (1999) to include stellar spiral arms, 
and showed that gas flows in models with stellar spiral 
arms match the observed I — v diagram better than models 
without stellar spiral arms. Fux (1999) performed three- 
dimensional TV-body + hydrodynamical simulations of the 
Milky Way, and gave a coherent interpretation of the main 
features standing out from observed I — v diagrams within 
the Galactic bar. In particular, he showed that the trac- 
ers of the gas associated with the Milky Way's dust lanes 
can be reliably identified and the 3 kpc-arm appears as 
a gaseous stream rather than a density wave (Fux 2001). 
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He also argued that the density centre of the stellar bar 
wanders around the centre of mass and the resulting gas 
flow is asymmetric and non-stationary. 

Advanced modeling of the Galactic disk provides fur- 
ther areas for comparisons between the theoretical models 
and observations, and, in addition, it can help us under- 
stand the true structure and dynamics of the Galactic 
disk. However, these previous numerical studies adopted 
a rather simple modeling of the ISM and the stellar 
disk. First, an isothermal equation of state (EOS) 
with a velocity dispersion of ~ 10 km s _1 (Wada et al. 
1994; Englmaier & Gerhard 1999; Fux 1999; Weiner & 
Sellwood 1999; Bissantz et al. 2003), or a phenomenolog- 
ical model of the ISM (Rodriguez-Fernandez & Combes 
2008) were used. The isothermal EOS would be a rele- 
vant approximation to investigate the global gas dynamics 
in galaxies, however the ISM is obviously not isothermal. 
It is important to include an energy equation with ap- 
propriate cooling and heating terms in order to compare 
the numerical results with the Hi and CO observations 
(e.g., Wada et al. 2000). The inhomogeneous nature of 
the ISM should be taken into account to investigate the 
gas dynamics on a local scale or in the central regions of 
galaxies (Wada & Koda 2001). In fact, the observed I — v 
diagrams (e.g. Dame et al. 2001) show many 'clumpy' 
sub-structures, which are not seen in the previous compu- 
tational I — v diagrams. 

Secondly, self-gravitational interactions of the ISM were 
often ignored (Englmaier & Gerhard 1999; Weiner & 
Sellwood 1999; Bissantz et al. 2003). The self-gravity of 
gas can play a significant role in high-density regions such 
as a central gaseous ring (Wada & Habe 1992; Wada et al. 
1994; Fukuda et al. 2000). Furthermore, the self-gravity 
in gas along with the thermal instabilities causes compli- 
cated, non-circular motions in the multi-phase ISM (Wada 
et al. 2002). A bar was modelled with a fixed potential 
undergoing a rigid rotation, and stellar spiral arms were 
not considered (Wada et al. 1994; Englmaier & Gerhard 
1999; Weiner & Sellwood 1999). In addition to the stellar 
bar, stellar spiral arms also result in non-circular motions 
of the gas (e.g., Fujimoto 1968; Roberts 1969; Shu et al. 
1973). Bissantz et al. (2003) included stellar spiral arms, 
but the spiral arms were assumed as rigid patterns. 

In order to understand the origin of the both large-scale 
structures and clumpy features in the observed I — v dia- 
grams, we present a self-consistent high-resolution simu- 
lation of a disk galaxy, which consists of a stellar disk and 
the multi-phase ISM in a static dark-matter halo. Star 
formation from cold, dense gas and energy feedback from 
type II supernovae (SNe), which have not been included 
in previous models, are taken into account. We have sim- 
ulated a model galaxy with a smaller mass than the Milky 
Way galaxy where the circular velocity at 8 kpc is 25 % 
less than that in the Milky Way galaxy. This enables us 
to ensure a high mass resolution for the ISM. Despite this 
difference, the / — v diagram from our simulations can sur- 
prisingly reproduces many features noted in observations, 
suggesting that our method has potential to qualitatively 
reproduce features in the Milky Way galaxy with a more 
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massive model. We are now preparing to run this 'Milky 
Way model' using ten times the number of particles (i.e. 
> 10 7 particles) used in the work presented here. 

In section 2, we describe our methodology and in sec- 
tion 3 we report the numerical results. We present the 
qualitative comparison for these structures of the Milky 
Way galaxy via I — v diagrams in section 4. In section 5, 
we give a brief summary. 

2. Numerical Method and Model Setup 

We used the N- body /hydro dynamic simulation code 
ASURA (Saitoh et al. 2008; Saitoh et al. 2009) to solve 
the Newtonian equation of motions and the equations 
of hydrodynamics using the standard smoothed particle 
hydrodynamics (SPH) methods (Lucy 1977; Gingold & 
Monaghan 1977; Monaghan 1992) : 

N nh 

Pi = ^2mjW(\xi-Xj\,h), (1) 

3 

+&-V$ DM (xi), (2) 




where m, p, p, u, v, x, and $dm are the mass, den- 
sity, pressure, specific internal energy, velocity, position 
of the gas, and the gravitational potential of the dark 
matter halo, respectively. We assume an ideal gas EOS 
p = (7 — l)pu, with 7 = 5/3. W(x,h) and h are the SPH 
smoothing kernel and the smoothing length, respectively, 
and h is allowed to vary both in space and time with the 
constraint that the typical number of neighbours for each 
particle is 7V n b = 32 ± 2. The artificial viscosity term 
(Monaghan 1997) and the correction term to avoid large 
entropy generation in pure shear flows (Balsara 1995) are 
used. Radiative cooling of the gas, A, was solved assum- 
ing an optically thin cooling function with solar metal- 
licity which covered a wide range of temperature, 20 K 
through 10 8 K (Wada & Norman 2001). Heating due to 
far-ultraviolet radiation (FUV) and energy feedback from 
SNe, T = Tpuv + Tsn, was also included. We assume a 
uniform FUV field with matches that observed in the so- 
lar neighborhood: 

Tfuv = 10- 24 eG n H [erg s^cm" 3 ] (4) 

(Wolfire et al. 1995), where e and Go are the heating ef- 
ficiency (e = 0.05) and incident FUV normalized to the 
solar neighborhood value (Go = 1), respectively. The self- 
gravity of stars and SPH particles, g, is calculated by 
the Tree with GRAPE method (Makino 1991). We here 
used a software emulator of GRAPE, Phantom-GRAPE 
(Nitadori et al. in preparation). 
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Models for the star formation and the SN feedback were 
the same as those in Paper I and Saitoh et al. (2008). We 
adopted the single stellar population approximation, with 
the Salpeter initial mass function (Salpeter 1955) and the 
mass range of 0.1 — 100 M . If an SPH particle satisfies 
the criteria (1) n H > 100 cm" 3 , (2) T < 100 K, and (3) 
V • v < 0, the SPH particle creates star particles following 
the Schmidt law (Schmidt 1959), with a local star for- 
mation efficiency, C* = 0.033, in a probabilistic manner 
(Tasker & Bryan 2006; Saitoh et al. 2008; Saitoh et al. 
2009; Tasker & Bryan 2008; Robertson & Kravtsov 2008). 
The local star formation efficiency C* is an unknown pa- 
rameter. Saitoh et al. (2008) however showed that global 
(galactic) star formation rate is not directly proportional 
to the local star formation efficiency, C*, but mainly con- 
trolled by gas mass in the high density regions, which 
is statistically related to the global evolution of the ISM 
(for details, see Section 5.2 in Saitoh et al. 2008). We 
treated Type II SNe feedback as injecting the thermal en- 
ergy from SNe to the neighbour SPH particles. We here 
assumed that stars with > 8 M experience Type II SNe 
and each SN release an energy with a canonical value of 
10 51 ergs. It was suggested that if the feedback energies 
are lower than the canonical value, it is hard to reproduce 
hot, diffuse halos around disk galaxies (Tasker & Bryan 
2006; Tasker & Bryan 2008). 

We follow the same method presented in Paper I: we 
first embed a pure iV-body stellar disk with an exponen- 
tial profile within a static dark-matter (DM) halo poten- 
tial whose density profile follows the Navarro- Frenk- White 
profile (Navarro et al. 1997). Using Hernquist's method 
(Hernquist 1993), we generate the initial conditions for 
the stellar disk which is then allowed to evolve for 2 Gyr, 
resulting in the formation of a stellar bar and multi-arms. 
At this point, we added a gaseous component to the stellar 
disk. This combined fixed dark-matter potential embed- 
ded stellar disk + gas disk is used as the initial condition 
for the simulation presented in this paper. The model pa- 
rameters of the dark matter halo, stellar disk, and gas disk 
are summarized in Table 1. The circular velocity curve is 
shown in Figure 3a. 

The total numbers of stars and gas particles are 3 x 10 6 
and 10 6 respectively, and particle masses are 11000 M 
and 3200 M . This gas particle mass allows us to resolve 
the gravitational fragmentation of dense clumps of gas 
down to ~ 10 5 M (Saitoh et al. 2008). In addition, a 
gravitational softening length of 10 pc in association with 
a typical smoothing length in dense regions of a few tens 
of pc means that structures larger than a few tens of pc 
are well resolved. 

3. Results 

3.1. Three-dimensional structures 

In Figure la, we present a face-on view of the stars at 
t = 1.24 Gyr where t = is the time when the gas compo- 
nent is added to the pre-evolved stellar disk. By this time, 
we see a stellar bar with a semi-major axis of ~ 3 — 4 kpc. 
In R = ~ 5-10 kpc, two- armed stellar spiral dominates, 
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Table 1. Model parameters for each mass component (dark 
halo, stellar disk, and gas disk) 



Component 


Parameter 


Value 


Dark Halo 


Mass 


6.3 x 10 11 M 




Radius 


122 kpc 




Concentration 


5.0 


Initial Stellar Disk 


Mass 


3.2 x 10 iU M 




Scale Length 


3.0 kpc 




Scale Height 


0.3 kpc 


Initial Gas Disk 


Mass 


3.2 x 10 y M 




Scale Length 


6.0 kpc 




Scale Height 


0.2 kpc 



and outside ~ 10 kpc multi-armed spirals appear (see also 
Figure 2 and 4 in Paper I). The Fourier amplitude of the 
dominant mode is ~ 0.1 — 0.2, which is consistent with ob- 
servational values of external disk galaxies (Rix & Zaritsky 
1995; Zibetti et al. 2009). 

Figure lb present a face-on view of the cold gas (T < 
100 K) overlaid on the stellar disk. Along the bar a typi- 
cal offset ridge structure of cold gas (dust lane) is formed 
and enhanced by the bar potential. There are many gas 
arms in the disk, but contrast to previous isothermal sim- 
ulations (e.g., Fux 1999), these arms are not smooth. 
There are also many substructures, such as clumps and 
filaments, in the inter- arm regions. Hence, it is not easy 
to trace any "single arm" from the central region to the 
outer disk, even for major spiral arms (see also Figure 7c). 
Typical length and width of the gas filaments are < 500 
pc and ~ 100 pc, respectively. These local structures, i.e., 
gas clumps and filaments, are originated from their self- 
gravity, whereas large-scale structures, i.e., gaseous spiral 
arms, are driven by the stellar bar and spiral arms. 

Figure lc shows the surface density of the stars on the 
(/,6)-plane. The observer is assumed to be at (i?o?0b) = 
(8 kpc, 25°), where Ro and 0b are the galacto-centric dis- 
tance and position angle relative to the major axis of the 
bar, respectively (positive sign means for the anti-galactic 
rotation). Since many previous studies suggest that the 
orientation of the major axis of the Galactic bar relative 
to us typically ranges from 15° to 35° (e.g., Binney et al. 
1997; Fux 1999; Englmaier & Gerhard 1999; Bissantz et al. 
2003; for a review, Merrifield 2004), we adopt a fiducial 
value of 0b = 25°. Due to our choice of 0b = 25° we observe 
an asymmetry in disk and bulge densities in longitude. 
This asymmetry agrees with near-infrared observational 
studies by COBE/DIRBE (Weiland et al. 1994; Dwek 
et al. 1995; Binney et al. 1997). 

Surface density contours of the stellar disk are shown 
in Figure 2. The principal axis of the surface density 
distribution for a given annulus of width AR = 0.5 kpc 
are indicated by thick lines. Within a radius R < 4 kpc, 
they well align in a same axis however, the outer annuli 
show a different distribution axis. Figure 3b shows pat- 
tern speeds of the principal axis. The pattern speeds for 
R < 4 kpc are more or less the same at ~ 27 km s _1 
kpc -1 , however further out the pattern speed declines. 
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Therefore, we determine that the semi-major axis of the 
stellar bar (Rb) is ~ 4 kpc with a face-on axis ratio of 
~ 0.5 and the bar rotates as a rigid-body with fib °° 27 
km s _1 kpc -1 . Comparing fib and the angular frequency 
of the galaxy (fi), we infer the location of the co-rotation 
radius (Rcr) be around R = 5 kpc, whose value is larger 
than R\). This is consistent with the theoretical stud- 
ies (e.g., Contopoulos 1980) where a self-consistent bar is 
required to have Rcn/Rb > 1, as well as observations of 
other barred galaxies (e.g., Aguerri et al. 2003). The inner 
and outer Lindblad resonances are located around R = 1 
kpc and 10 kpc, respectively. We note that the values of 
these bar parameters do not significantly change during 
the simulation (t = 2 Gyr). The inner Lindblad resonance 
(ILR) locates at ~ 1 kpc in this model, which is quite 
different from some previous models (e.g., Binney et al. 
1991; Bissantz et al. 2003) that place the ILR at - 100 
pc. This difference may be caused by differences of the 
mass model in the central region between our model and 
the Milky Way (see section 4). 



TIME = 1 .300 Gyr 




-6 -4 -2 2 4 6 
X [kpc] 

Fig. 2. Surface density contours of the stellar disk at t = 1.3 
Gyr are spaced by a constant interval of 0.2 dex. Directions of 
the principal axis from R = 0.5 kpc to 5.0 kpc every 0.5 kpc are 
indicated by thick lines. The galactic rotation is clock-wise. 




R [kpc] 

Fig. 3. (a) Circular velocities of the galaxy model. The cir- 
cular velocities (solid curve) are derived from the azimuthally 
averaged potential. The contribution from dark matter halo 
is shown by a dashed curve, (b) Angular frequencies as a 
function of radius. Rotation angular (Q) and epicyclic (/c) fre- 
quencies are derived from the azimuthally averaged potential. 
The dashed curves indicate Q — l/4«. Thick horizontal lines 
are pattern speeds of the principal axes presented in Figure 2. 



3.2. The I - v features 

Figure 4 shows a synthetic I — v diagram derived from 
the cold gas distribution given in Figure la. The diagram 
shows not only the peculiar features but also the terminal 
velocity tangent points observed in the H I and CO I — v di- 
agrams. Furthermore, our synthetic l — v diagram shows a 
clumpy morphology whose features posses a large velocity 
width (typically ~ 10 — 20 km s _1 ), in agreement with CO 
I — v diagrams (Dame et al. 2001). These clumpy features 
have not been reproduced in the previous studies (Wada 
et al. 1994; Englmaier & Gerhard 1999; Bissantz et al. 
2003; Fux 1999; Rodriguez-Fernandez & Combes 2008). 





I [deg] 

Fig. 1. (a): V-band image of the stellar disk (orange) at t = 1.24 Gyr. A major axis of the stellar bar inclines toward 25° from a 
y-axis. The galactic rotation is clock- wise, (b): Overlaid of the cold gases (green, T < 100 K) on the panel (a), (c): Distribution 
of stellar surface density on the l — b plane. The observer is located at the red point (0 kpc, — 8 kpc) in the panels (a) and (b). 



We also show the I — v diagram from a simulation with an 
isothermal EOS 2 in Figure 5a. Here this model includes 
the self-gravity of the ISM. Unlike the I — v diagram of the 
run with the multi-phase ISM, this I — v diagram shows a 
smooth rather than clumpy morphology. This result in- 
dicates that the appropriate treatment of the gas physics 
in ISM, especially cooling down to T <C 10 4 K, is essential 
to understand clumpy morphology formed in the observed 

1 — v diagrams. 

The features in I — v diagrams depend on both the 
location of the observer within the galactic plane and 
the time. To illustrate this, we show I — v diagrams for 
GRo,0b) = (8 kpc, 45°) and (6 kpc, 25°) in Figure 5b and 
5c, respectively. If we observe the nuclear region of the 

2 We replaced the gas in the multi-phase run into the isothermal 
gas with T = 10 4 K at t = 1.0 Gyr, and then let it evolve for 
240 Myr. Thus the stars distribute almost same as in the multi- 
phase run. 



simulation from 0b = 45°, the velocity of the core is lower 
than that for the case with 0b = 25°. The diagonal fea- 
tures, such as so-called '3-kpc arm' (labeled by E in Figure 
6) seen in the CO I — v diagrams, also become shallower. 
The position of the 'Carina arm' (labeled by F in Figure 
6) around I ~ —50° — 60° and v ~ 30 km s _1 also shifts 
to larger I. Changing the radial position of the observer 
to R = 6 kpc alters the I — v diagrams such that we see a 
larger number of high- velocity components in the regions 
of I < 0,v > or I > Q,v < 0. The l — v diagram at t = 1.34 
Gyr, which is 100 Myr after the snapshot used for produc- 
ing Figure 4, is shown in Figure 5d. One can notice that 
the non-axisymmetric features change on this time-scale. 
This reflects the fact that the stellar and gaseous spiral 
arms are not stationary, but time-dependent phenomena. 
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longitude [deg.] 

Fig. 4. Synthetic I — v diagram in the range \b\ < 1° derived from the cold gas (T < 100K) distribution given in Figure la (t = 1.24 
Gyr). The observer is located at the red mark point indicated in Figure 1, and has pure circular motion. The contribution of each 
SPH particle to the synthetic I — v diagrams is weighted by its inverse squared distance relative to the observer in order to mimic 
the flux decline of point source. The gases with d < 500 pc are not displayed. 



4. A comparison to real structures in the Milky 
Way 

4-1. Large- Scale Structures 

A qualitative comparison between the numerical results 
and observations of the Milky Way is given in Figure 6. 
We identify and trace schematically (solid lines) such pe- 
culiar features as the Perseus (A), Outer (B), connecting 
(C), 135-kms- 1 (D), 3-kpc (E), and Carina (F) arms. The 
terminal velocity tangent points such as the Sagittarius 
(a), Scutum (b), Norma (c), and Centaurus (d) tangent 
points, are also shown schematically by squares. Although 
the Perseus arm, Sagittarius tangent point, and central 
molecular zone (CMZ) are unclear in our model, other 
structures are clearly evident. 

In order to see which features in the I — v diagram cor- 
respond to structures in the real space, we marked several 
major features with the same color in the face-on map 
(Figure 7a, b). This comparison gives us insight to in- 
fer the real morphology of spiral arms in the Milky Way. 
For example, a prominent feature known as the "molecu- 
lar ring" observed in CO l — v diagrams is not necessarily 
correspond to a "ring" in real space, but be a part of 
nearby spiral arm (such as the the one colored by light 
blue). Nakanishi & Sofue (2006) suggest that the molecu- 
lar ring results from a combination of the inner part of the 
Sagittarius-Carina arm and the Scutum-Crux arm. The 
"3-kpc arm"(E) corresponds to an inner spiral arm (the 
one colored by green). It has been proposed that this arm 
arises from a lateral arm surrounding the bar (Figure 16 
in Fux 1999), or a small arm starting from the bar-end 
(Figures 13 and 15 in Bissantz et al. 2003). Following the 
interpretation by Bissantz et al. (2003), the 3-kpc arm 



would correspond to the gaseous arm colored in blue in 
Figures 7a, b. Therefore, our result supports to the in- 
terpretation by Fux (1999). The near part of the offset 
ridge is the one colored red, which could be due to a part 
of the "connecting arm"(C). Therefore, the clumpy na- 
ture in this arm means that the offset ridge of the Milky 
Way galaxy can consist of some gas clumps. Contrary to 
Fux (1999), "135-km s _1 arm" (D) is a part of the far 
side of the bar end in our model. We note, however, that 
above arguments on this comparison between correspond- 
ing structures in our model and observations is qualitative, 
because our current "barred" galaxy is somewhat smaller 
than the Milky Way resulting in a lower rotational velocity 
at 8 kpc of 163 km s _1 compared to ~ 220 km s -1 . 

The cold gas in the galactic center (< 200 pc, the one 
colored by purple) corresponds to the feature known as 
CMZ in the observed CO I — v diagram. However, it is 
not clearly seen in our I — v diagram model in comparison 
to that in the Milky Way galaxy. On the mass modeling, 
this may be attributable to differences in the bar parame- 
ters between our model and the Galactic bar, but also to 
absence of a central bulge. It is known that bar parame- 
ters, such as the pattern speed and axis ratio, can affect a 
gas inflow rate and shape of dust lanes, thereby influence 
formation of the CMZ (Athanassoula 1992; Englmaier & 
Gerhard 1997; Patsis & Athanassoula 2000). According to 
these previous studies, the Galactic bar may be rounder 
or more concentrated than the bar obtained in our simu- 
lation. Another discrepancy between our model and the 
Milky Way is that stellar number counts show evidence 
for the existence of an inner bar in the Milky Way, which 
is much smaller than the outer bar (or triaxial bulge) with 
a semi- major axis of 3.5 kpc (Alard 2001; Nishiyama et al. 
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Fig. 5. Same as Figure 4, but for (a) the isothermal gas with 
T = 10 4 K, (b) b = 45°, (c) Rq = 6kpc, and (d) t = 1.34 Gyr. 



2005) 3 . However, our model does not have nested stellar 
bars. Rodriguez-Fernandez & Combes (2008) suggest that 
the parallelogram shape of the I — v diagram of the CMZ 
is due to the influence of this inner bar (also known as nu- 
clear bar or secondary bar). The central bulge, which is 
not included in the initial condition of the present model, 
may also affect the dynamics of the central region. More 
detailed comparison between the observed kinematics of 
the gas and models in the central several degrees requires 
high numerical resolution in the central part in numerical 
models. Effects of the central bulge and bar properties 
on the kinematics will be discussed in the high resolution 
models elsewhere. 

In terms of realistic modelling heating physics of the 
ISM, we need to model the spatial dependence of the 
incident FUV, which has been assumed to be constant 
(Go = 1.0) everywhere in this paper (see equation 4). 
Gerritsen & Icke (1997) calculated the incident FUV by 
summing the FUV flux from all stars according to their 
ages. The resultant radiation field declines outside the 
stellar disk with radius as 1/R 2 . This means that the in- 
cident FUV is stronger in the central region than that in 
the solar neighborhood, implying that if we were to con- 
sider the spatial dependence of the incident FUV in our 
model, for example Go oc 1/R 2 , then the consumption of 
gas by star formation in the region of the CMZ might be 
suppressed, yielding a stronger CMZ present in our I — v 
diagram. 

4-2. Clumpy Morphology 

In Figures 7c and 7d, we marked gas clumps (> 10 5 M ) 
with the same color between the face-on map and I — v di- 
agram. We here defined the gas clumps by using a Friend- 
of-Friend (FOF) method (Davis et al. 1985) with a linking 
length of 30 pc for the cold gas. The large-scale pecu- 
liar features, such as the "Connecting arm"(C), "3-kpc 
arm"(E) and "135-km s _1 arm"(D) in I — v diagram are 
actually ensembles of dense clouds and filaments. 

It is known that the Bania's "Clump 1" and "Clump 2" 
are placed around the negative longitude end of the 135 
km s _1 arm and at the positive longitude side of the CMZ 
in observed CO I — v diagram (Bania 1977; Bania et al. 
1986; Stark & Bania 1986). We associate these clumps 
with clumps "cl" and "c2" in Figures 7c and 7d, respec- 
tively. Mass of the clump "cl" is ~ 10 5 M , the clump 
"c2" is ~ 5 x 10 5 M . From time evolution of these 
clumps, we found that the clumps "c2" eventually col- 
lapsed into the galactic center. Other gas clumps within 
a bar region tend to spiral into the galactic center with 
losing their angular momenta around their orbital apocen- 
ters via collisions between them. Fux (1999) inferred that 
the Bania's clump 2 and another vertical feature near the 



There is a discussion on the existence of the inner bar in the 
Milky Way. At the central region, the magnitude differences 
in the star counts on both sides are less than 0.1 mag, so ef- 
fects from asymmetric dust absorption will cause quite some 
uncertainly in the interpretation of the star counts. However, 
observations of external barred galaxies shows that many barred 
galaxies have inner bars (e.g., Erwin 2004). 
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clump could be correspond to the gas clumps which are 
just going to cross the near-side offset axis (i.e. connect- 
ing arm). Contrary to this picture, we here suggest that 
vertical features in observed CO I — v diagrams near the 
CMZ are streams of gas clumps spiraling into the galactic 
center. As a consequence, our simulation suggests that gas 
is stochastically supplied into the galactic central region 
by dense gas clumps. 

4-3. The Non-Stationarity of l — v Features 

We find that the observed I — v diagram is reproduced 
only at a specific time (t ~ 1.24 Gyr) and their features 
change on a time-scale of ~ 100 Myr, suggesting that 
the observed I — v features are transient features. The 
non-stationarity has been also found by Fux (1999) and 
Bissantz et al. (2003) in their numerical models. Fux 
(1999) and Bissantz et al. (2003) attributed the non- 
stationarity to wandering of a stellar bar around the cen- 
ter of the mass, and the difference of the pattern speeds 
between the stellar bar and spiral arms, respectively. In 
fact, our model shows that the stellar spiral arms rotate 
slowly than the bar (Figure 3b), suggesting that decou- 
pling between the stellar bar and spirals may cause the 
non-stationarity. 

Recently, Baba et al. (2009) analyzed kinematics of star 
forming regions in the same numerical model discussed in 
this paper with observed proper motions of maser sources 
in the MW, and suggested that the Galactic stellar spiral 
arms should be transient, recurrently formed structures 
rather than the "stationary" density waves proposed in 
Lin & Shu (1964) (for a review, see Bertin & Lin 1996). 
This transient nature of stellar spirals is also supported by 
previous TV-body simulations of stellar disks without the 
ISM (Sellwood & Carlberg 1984; Sellwood 2000; Sellwood 
2010; Fujii et al. 2010). We infer that dynamic nature of 
stellar spiral arms themselves could contribute the non- 
stationarity of the I — v diagram. We will quantitatively 
investigate a driving mechanism of the non-stationarity in 
our following paper. 

5. Conclusions 

By using high-resolution, 7V-body+hydrodynamical 
simulation in which the multi-phase ISM, star-formation, 
and SN feedback were self-consistently taken into account, 
we qualitatively reproduced not only large-scale structures 
of the Hi and CO I — v diagrams such as the terminal ve- 
locity tangent points and the coherent features, but also 
clumpy structures. Previous studies with numerical simu- 
lations on I — v diagrams did not reproduce these clumpy 
structures. When we adopt a model galaxy whose ve- 
locity is similar to the Milky Way galaxy with the same 
numerical method, we can advance our argument more 
qualitatively. We will show the results in the near future. 

The authors are grateful to the anonymous referee for 
his/her valuable comments. We would like to thank 
James Binney for his constructive comments for our 
previous paper, by which led us to analyze our re- 
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Fig. 7. (Left panels) Correspondence between the gaseous spiral arms in the x — y plane and the I — v features for the cold gas in 
the inner galaxy at t = 1.24 Gyr. (a): Spacial distribution on the x — y plane, (b): l — v diagram. Each SPH particle is plotted as a 
dot in spite of its distance. The observer is located at (Ro,4>b) — (8 kpc, 25°) presented by the red mark point in the top panel, and 
has a pure circular motion. (Right panels) Same as the left panels, but colored for gas clumps. Solid curves labeled by "C", "D", 
and "E" correspond to them in Figure 6. Open circles (cl and c2) are examples of gas clumps around their apocenters of orbits. 



